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In this paper we study the kinetics of diffusion-hmited, pseudo- first-order A + B ^ B reactions in 
situations in which the particles' intrinsic reactivities are not constant but vary randomly in time. That is, 
we suppose that the particles are bearing "gates" which fluctuate in time, randomly and independently of 
each other, between two states - an active state, when the reaction may take place between an A and a B 
particles appearing at a close contact, and a blocked state, when the reaction is completly inhibited. We 
focus here on two customary limiting cases of pseudo-first-order reactions - the so-called target annihilation 
and the Rosenstock trapping model, and consider four different particular models, such that the A particle 
can be either mobile or immobile, gated or ungated, as well as ungated or gated B particles can be fixed 
at random positions or move randomly. All models are formulated on a d-dimensional regular lattice and 
we suppose that the mobile species perform independent, homogeneous, discrete-time lattice random walks. 
The model involving a single, immobile, ungated target A and a concentration of mobile, gated B particles 
is solved exactly. For the remaining three models we determine exactly, in form of rigorous lower and upper 
bounds showing the same A^-dependence, the large- asymptotical behavior of the probability that the A 
particle survives until the A-th step. We also realize that for all four models studied here the A particle 
survival probalibity can be interpreted as the moment generating function of some functionals of random 
walk trajectories, such as, e.g., the number of self-intersections, the number of sites visited exactly a given 
number of times, "residence time" on a random array of lattice sites and etc. Our results thus apply to the 
asymptotical behavior of the corresponding generating functions which has not been known as yet. 
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I. INTRODUCTION 



Many of naturally occuring chemical reactions, or reactions used in various technolog- 
ical and material processing operations involve molecules with a rather complex internal 
structure. For such reactions the geometrical complexity of the molecules appears to be 
a significant rate-controlling factor, additionally to the transport limitations and the con- 
straints imposed by the elementary reaction act; that is, chemically active groups of complex 
molecules involved may be effectively screened by its inactive parts, which thus impedes the 
access of the reactive species and inhibits the reaction. For instance, geometrical restric- 
tions are crucial for ligand binding to proteins, such as, e.g. myoglobin or hemoglobin [|l|. 
Here, in the static X-ray structure of myoglobin there is no hole for the ligand to enter 
and it is believed that binding of the ligand happens when the side chains blocking the 
entrance swing out in course of their thermal motion (Figs.l and 2). Similarly, intercalation 
of drugs by DNA may be controlled by breathing motions that involve the unstacking of 
adjacent pairs of bases. In some cases, ligands themselves can possess a complicated internal 
structure, as exemplified, for instance, by the peptides, such that their reactivity may be 
influenced significantly by conformational changes. Last but not least, geometrical restric- 
tions do manifest themselves in such contexts, as, e.g., certain chemical reactions occuring 
within biomembranes 0, incoherent exciton trapping by substitutional traps on aromatic 
vinyl polymers molecular transport inside proteins 0-0, as well as in some medical 
therapies f^. Clearly, understanding of the impact of the geometrical limitations on the 
reaction kinetics constitutes an important challenge for the theoretical analysis. 

A physically plausible approach to account for the influence of the geometrical restrictions 
on the reactivity of complex molecules is to assume that the reaction in question, say, a 
generic pseudo-flrst -order reaction of the form 

A + B-^B, (1) 

is modulated by the side reactions of the form A ^ A* ot B ^ B*, where A (B) stands 
for an active, while A* (B*) - for an inactive, blocked state, in which case the reaction in 
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Eq.(|l]) is complete inhibited. In other words, one says that one or both species involved 
in the reaction are being gated, the gates changing in time their states according to some 
prescribed rules. 

Kinetics of gated diffusion-limited reactions has been studied analytically for nearly two 
decades. Following the seminal work of McCammon and Northrup who have analysed 
a simple case of a two-state gating described by arbitrary deterministic function of time, 
several important advancements have been made. In particular, classical Smoluchowski ap- 
proach |jTo| has been generalized in Refs. |ini-[l3|| to describe kinetics of stochastically-gated 



(SG) pseudo-first-order reactions in Eq.([l|) for both cases when the gate is imposed on the 
traps B or on the A particles. It has been realized that no symmetry exists between these 
two situations; as a matter of fact, the kinetic behavior appears to be rather different de- 



pending on which of the species precisely is being gated. Furthermore, Spouge et al |T4 
studied the SG reactions with more general mechanisms, including a non-Markovian case, 
while Berezhkovski with collaborators discussed the impact of the many-particle effects on 



the SG reactions kinetics |TB|,|T3|,|TC[ . However, the available theoretical analysis is either 
based on uncontrollable assumptions, akin to the Smoluchowski-type approaches (see, e.g. 
|T^,|TB| and references therein), or employs exact formalisms, which do not always allow for 
explicit calculation of the corresponding decay patterns and become tractable only when 
some simplifying assumptions are made. Consequently, except for the relatively simple 
model involving an immobile ungated target A and a concentration of mobile SG B parti- 
cles (see, e.g. Refs. [|lT]-|T^), temporal evolution of the SG reactions remains incompletely 



understood. This is, of course, not at all surprising since diffusion-limited reactions with 
stochastic reactivity clearly pose more complex technical problems to the theoretical anal- 
ysis than their ungated counterparts, which themselves are not solvable exactly and often 
exhibit spectacular, essentially many-particle behavior. 

In this paper we study in detail the kinetics of the SG pseudo-first-order reactions in 
Eq.(|T]), involving a single A particle and a concentration of B particles. We focus here on 
two limiting cases - the so-called target annihilation ||l9|^3i and the Rosenstock trapping 



model pj], and consider four different models such that an A particle can be either mobile 
or immobile, gated or ungated, as well as ungated or gated B particles can be fixed at 
random positions or move randomly (Fig. 3). For computational convenience, all models are 
formulated on a d-dimensional regular lattice and we suppose that mobile species perform 
independent, homogeneous, discrete-time lattice random walks. Further on, in regard to the 
reactivity fluctuations we restrict ourselves to the two-state Poisson gating model of Ref. 
in which each gate is supposed to be in either of two states - an active and a blocked 



one, and to update its state at each tick of the clock. The updating process is assumed to 
proceed completely at random, without memory in time and without correlations with the 
gates imposed on the other particles. 

For such a gating model we analyse the time evolution of the A particle survival proba- 
bility P/V; i-e. the probability that a single mobile or immobile, gated or ungated A particle 
is not annihilated up to the A^-th time step by a concentration of immobile or mobile, un- 
gated or gated B particles. For the model involving a single immobile, ungated target A 
and a concentration of SG mobile B particles the complete temporal evolution of P/v is 
calculated exactly. For the remaining three models we determine exactly the form of the 
large- asymptotical decay of by deriving rigorous lower and upper bounds, which both 
show the same A^-dependence but slightly differ in prefactors. 

An interesting by-product of our analysis, which appears in especially lucid fashion within 
the discrete-space description, is that for any particular model the A particle survival prob- 
ability P/v can be interpreted as a moment generating function of certain functionals, which 
mirror the internal geometry of random walks trajectories - an issue which has recently 



gained a renewed attention in view of many important applications [p5H28[|. As a matter 
of fact, it has been known for a long time that for ungated trapping reactions taking place 
on a d-dimensional lattice the survival probability Pjv can be thought off as the moment 
generating function of the number iS(^{riv}) of distinct sites visited by a given trajectory 
{rjv} of an A^-step lattice random walk as well, for Brownian motion in a (i-dimensional 
continuum Pjv can be interpreted as the generating function of the volume of the so-called 



Wiener sausage, i.e. the volume swept by a spherical particle during time t ||3^. We realize 
that for the SG pseudo-first-order reactions some other characteristic functionals of random 
walk trajectories come into play, which probe some interesting aspects of the geometrical 
structure of a single or of a collection of lattice random walks. Depending on the particular 
model P/v appears then as the generating function of such characteristic functionals of ran- 
dom walk trajectories as, e.g., the number of self- intersections, the number of sites visited 
exactly a given number of times, "residence time" on a random array of lattice sites, as well 
as some others. Consequently, our results also apply to the asymptotical behavior of the 
corresponding generating functions, which has not been known as yet in many cases. 

The paper is structured as follows. In section II we discuss two customary limiting cases 
of the psedo-first-order reactions in Eq.(l) - the so-called target annihilation model and the 
Rosenstock trapping model, and present a brief outline of different results concerning their 
kinetic behavior in the ungated case. In section III we consider a gated target problem, 
focusing on the survival probability P/v of an immobile ungated target in the presence 
of randomly moving, stochastically-gated B particles. We derive here an exact temporal 
evolution of the probability that none of mobile particles hits the target within an A^-step 
walk. Further on, in section IV we study the survival of an immobile gated target in the 
presence of a concentration of mobile particles. The corresponding survival probability is 
found here in form of rigorous lower and upper bounds, which both display the same A^- 
dependence and thus determine the temporal evolution of P/v exactly. Next, in sections 
V and VI we consider two models of stochastically-gated trapping reactions; in the first 
case these are the immobile traps B who are supposed to be stochastically gated, while the 
mobile A particle is assumed to be always in the reactive state (ungated), and in the second 
case the mobile A particle is assumed to bear a stochastic gate, while the immobile traps B 
are being considered as perfect, non- fluctuating traps. For both cases the A particle survival 
probability is determined exactly, in form of rigorous lower and upper bounds showing the 
same A^-dependence. Finally, we conclude in section VII with a brief summary of our results. 
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II. A REMINDER ON THE KINETICS OF UNGATED PSEUDO-FIRST-ORDER 
REACTIONS: TARGET ANNIHILATION AND THE ROSENSTOCK TRAPPING 

MODELS. 



To fix the ideas, we present first a brief summary of results concerning tlie kinetics of un- 
gated, diffusion-limited pseudo-first-order reactions in Eq.(|l]) taking place on (i-dimensional 
regular lattices. We will focus here and in the remainder of the paper on two particular 
cases - the so-called target annihilation model, involving a single, immobile A particle and a 
concentration of randomly moving B particles, and the Rosenstock trapping model involving 
a single mobile A particle performing lattice random walk in the presence of a concentration 
of immobile, randomly placed B particles - the traps. In subsequent sections of our work we 
will study how the kinetics in these two models will be modified if one imposes a stochastic 
gate on either of two species; consequently, the results of this section will serve us in what 
follows as an important point of reference. 



A. Target annihilation. 



We start with the target annihilation model which allows for an exact solution |T^-^. 
Consider an immobile A particle, located at the lattice origin, and B particles, which are 
initially all placed at random positions on a d-dimensional regular lattice and after that 
are allowed to move by performing independent, homogeneous, discrete-time random walks. 
As soon as any of the Bs appears at the lattice origin, the A particle gets instantaneously 
annihilated. Thus, following the earlier introduced terminology, the A particle will be called 
here as the "target", while the B particles will be referred to as the "scavengers" [^]. The 



property whose temporal evolution we wish to study is the probability Pn = exp[— cQ 



(tar)-, 
N J 



that the target particle "survives" until time N; here and henceforth c will denote the number 
density of the B particles, while Q^"^'' , in view of the pseudo-first order of the reaction in 
Eq.(P, can be interpreted as the integral effective reaction rate. 
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Now we define the model more precisely. We first suppose that the lattice is of a finite 
extent and contains M sites, whereas the number of the scavengers is also fixed and is equal 
to K. In what follows we will turn to the limit M, K —>■ oo, while keeping the ratio K/M 
fixed, K/M = c. 

Further on, let S^'^-' denote the position at which the k-th scavenger appears on the n-th 
step (n = 0, 1, . . . , N) for a given realization of its trajectory. Then, we construct a function 
\1/Ar of the form 

N K 



N 



lim exp 

N K 



(fc) 



n=0 fe=l 



n=Ofc=l 



(2) 



where T (S^'^) j is the indicator function. 



X(X) 



(3) 



1, X = 0, 
0, X^O, 

which shows whether the A;-th scavenger is at the lattice origin at the n-th step, or else- 
where. 

Note, that the function serves as the indicator of the reaction event; it is equal to one 
if within an A^-step walk none of the K scavengers has visited the origin, (i.e. the target), 
and turns to zero when within the "time" interval [0, A^] at least one of the scavengers has 
at least once visited the origin. 

The property one is generally interested to compute is not, however, the realization- 
dependent function "^n, but rather its averaged value Pn =< "^n >, the average being 
taken over different realizations of scavengers' random walks and their initial positions. 
Below we will briefly outline an exact computation of < \l'jv > (see also Refs. |T9|-p^j). 

Noticing first that all S^'^^ with different values of k are independent of each other, we 
have that the target survival probabilitity can be written in the following factorized form 

(4) 



Pn = (lim^^ooTj E ^So { exp [ - /3 ^ J (S„) ] } 

^'^ So n=0 



where the summation extends over all sites of a d- dimensional lattice, while the symbol 
Esp| ... I denotes the expectation on the set of different random walk trajectories starting 
at the site Sq. Note next that 

TV 

Eso{lim^^ooexp[-/3 J] X (S„) ] } = Prob(S„ ^ for any n G [0, iV]|So) , (5) 

n=0 

i.e. is equal to the probability that a random walker starting at the site Sq does not visit the 
origin within first steps. Turning now to the infinite-space limit, i.e. letting K, M ^ oo 
while keeping their ratio fixed, K/M = c, we find then that Eq.(^ attains the form 

\ ' II _ Prnb(S„ ^ n fnr anv n. 10 /VIIS^l 



Pn = limK,M^oo\K/M=c{} - H ~ Prob(S„ ^ for any n e [0, A^]|So] 

So 

c^Prob^(o|So)], (6) 



exp 







where ProbAr(^0|So) stands for the probability that a first passage from the site So to the 
origin did actually occur within the first steps. 

The probability Prober (^0|So) is known exactly (see, e.g. |^); being summed over all 
possible starting points it defines another important characteristic property of random walks 
- the expectation of the number 5(^{SAr}^ of distinct sites visited by a given trajectory {Sat} 
of an A^-step walk starting at the origin, i.e., 

5:Prob;v(o|So) = Eo{5({S^})} (7) 

So 

The expected number of distinct sites visited by an A^-step walk shows different asymptotical 
behavior depending on the dimensionality d and on the type of the lattice pSf . In particular, 
for (i- dimensional Polya random walks one has 



,.1 E„{5({S„})}^(?^)'%0(-^) 



A^^ 

d = 3 E„{s{lS,})}^-p^^ + 0(s/N), (8) 

where C2 = 4/3 -\/3, 1 and 2/-\/3 for hexagonal, square and triangular two-dimensional 
lattices, respectively. The constant P(0|0; 1~) determines the probability R of eventual 



return to the origin, P(0|0; 1~) = [1 — Rj , and is defined as the hmit ,^ ^ 1^ of the 
generating function P(0|0;^) = J2n=o -^^{^10)^^ , PAr(0|0) being the probabihty of finding a 
random walker at the origin at the iV-th step, provided that the walk has started at the origin. 
The exact values of P(0|0; 1") are 12T%1)/{2^/'^tt^), VQ r(-^) r(^) r(^) r(ii)/(284 vr^), 
r^(i)/(4 n^) and r®(i)/(2^^/^7r'^) for diamond, simple cubic, body-centered cubic and face- 
centered cubic lattices, respectively [p5|] . 

Therefore, for the target annihilation model the decay law can be computed exactly and 
the integral effective reaction rate Q^^^^^ equals simply the expected number of distinct sites 
visited by an A^-step walk [p!9|-p2|. Note also that the decay form appears to be essentially 
dependent on the dimensionality of the embedding lattice; it is characterized by a stretched- 
exponential dependence for low dimensional lattices, on which the Polya walks are recurrent 
[R = 1,P(0|0;1~) = oo), and shows a purely exponential behavior for lattices of spatial 
dimension d > 2, where R < 1 and P(0|0; 1~) is well defined. 



B. Rosenstock trapping model. 

We turn next to the so-called Rosenstock trapping problem, in which one focuses on 
the fate of a single A particle performing a random walk over the lattice in the presence of 
immobile, perfect, randomly placed traps B. 

We start assuming again that the lattice is finite and contains M sites. The K traps B 
are placed completely at random and their positions are determined by the lattice-vectors 
S^'^-', k = 1,2, . . . , K. Denoting the lattice position of the A particle at the n-th step as r„, 
we can now write down the indicator function of the reaction event as follows 

N K 

= lim/3^oo exp[-P 2:(r„ - S^^))] , (9) 

n=Ok=l 

where ^{p^ is the indicator function defined in Eq.(|^). The function in Eq.(P) is equal to 
one for such A^-step trajectories which avoid passing through any of the sites S*-*^-* and turns 
to zero for those trajectories which visit at least once at least one of these sites. 
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Now, we pass to averaging the function in Eq.(|^) over the traps' placement. Since all 
g(fc) g^j^g mutually-independent, one can wright down such an the average in the factorized 
form 

K 1 N 



1 

< ^jv > = Eo{ n (t7 E lim/3^oo exp 
k=i ^'^ s{fe) 



/?E^(r.-S('=))])} 



n=0 



(10) 



Next, in the limit K, M ^ oo one has that 



/3^oo exp 



/Ji:x(r„-S)])^-} 



n=0 



Eo{ exp 
Eo{ exp 



N 

^ (l - lim^^oo exp - (3 x(r„ - S) ) } 

n=0 



N 

E(l-x(El(r„-S)))l}, 

n=0 



(11) 



where 



N 

E(i-^(E^(rn-s)))^5({M), 



(12) 



n=0 



is the number of distinct sites visited by a given trajectory {rjv}. Consequently, the A 
particle survival probability obeys 



< \E'Ar >= Eo| exp -cS[{ym]) | = exp -cQ 



(13) 



and hence, appears here as the moment generating function of the number of distinct sites 
visited by an A^-step random walk. 

Therefore, the major difference between the target annihilation model and the trapping 
model is exactly that in the former the survival probability is the exponential of the expected 
number of distinct sites, while in the latter case it involves a fairly more complex property 
- its moment generating function. In consequence, the trapping problem turns out to be 
essentially more difficult than the target one and hence, shows a richer behavior. 

To display the time evolution of Pn in the trapping problem, we will first outline the 
predictions of a certain heuristic approach - the so-called Rosenstock approximation 



and then write down the results of a rigorous analysis by Donsker and Varadhan, which 
concerns the large- asymptotical behavior ^ . 
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1. Rosenstock approximation. 



This approximation has been first apphed by Rosenstock in his studies of lumines- 
cence quenching kinetics and amounts, in essence, to the mere replacement of the average 
of an exponential of the number of sites visited by an exponential of the expected number. 



< "^N >~ exp 



cEo{5({r;v})}] (14) 

As one may readily notice, this heuristic procedure yields exactly the behavior found for the 
target annihilation problem and consequently, within the framework of this approximation 
one finds that Q^^^ = Q^l^^\ 

As a matter of fact, numerical simulations demonstrate that a rather crude and un- 
controllable approximation in Eg . (11^) provides quite a fair description of the decay for the 
trapping problem for intermediate values of [0, until at very large A^ some deviations 
emerge. The reason why the Rosenstock approximation works at intermediate A^ can be ap- 
parently explained as follows: As a matter of fact, the Rosenstock approximation represents 
a rigorous lower bound on Pjv? since replacement of the average of an exponential of the 
number of distinct sites visited by an exponential of the expected number is tantamount to 
the apphcation of the Jensen inequality, < exp[— c Qn]>^ exp[— c < Qn >]• On the other 
hand, this inequality can be rewritten as 



<\I/Ar>=<exp ln(\l/Ar) >> exp < ln(\E'Ar) > 



(15) 



which signifies that in such an approach the average of the idicator function is approximated 
by the exponential of the averaged logarithm of this function. Since logarithm is a slowly 
varying function, it is generally believed that its average value is supported by typical 
realizations of random walk trajectories which are representative at moderate values of A^. 



2. Fluctuation-induced large-N tails of the survival probability. 



In the large- A^ limit, however, kinetics of the trapping reactions proceeds somewhat 
slower than that predicted by Eq.(|lJ). This happens, namely, due to some fluctuation 
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effects, which a mean-field- type approximation in Eq.([T^ can not capture. It has been first 
predicted in Ref. and subsequently proven by Donsker and Varadhan that for 
arbitrary d the decay follows 

Pjv ~ exp[-arfc2/('^+2)Ar'^/('^+2)], N oo, (16) 
where ad is a constant, dependent on the lattice dimensionality pO |. 



The physical origin of such an anomalous decay law can be illustrated by the following 
heuristic derivation. Consider first the function in Eq.(|l^) and suppose that for some given 
realization of the A^-step A particle trajectory {yn} the maximal excursion from the origin 
is equal to Rmax- Consequently, for this realization the number iS(^{rAr}^ of distinct sites 
visited by this realization of random walk trajectory can be majorized as iS(^{r7v}) < Rmaxi 
and the overall decay function can be bounded from below by 

Pn> J d'^Rmax exp -cR'l^^^ Prob (max { rjv } = i?„„^), (17) 

where Prob(maa;{rAr} = Rmax) is the probability that for a (i-dimensional, A^-step random 
walk the maximal displacement from the starting point is exactly equal to Rmax- For suffi- 



''max 



where 



ciently large iV, the leading behavior of this probability follows ~ exp — 7dA^/ R^ 
■jd is a dimension- dependent constant. Substituting the latter form to Eq. ([l7|) and noticing 
that the integrand is a bell-shaped function, we thus perform the integration in terms of 
the saddle-point method. In doing so one finds that the value of Rmax which provides the 
maximum to the integrand is given by Rmax = {2'jdN/cdY^^'^'^'^\ Consequently, the overall 
decay function obeys 

Pn > exp [ - const c^'^''^^^ iV^/C^+s)] ^ ^^g^ 

which bound displays exactly the same A^-dependence as the result in Eq.(0). 

Another illustrative derivation can be performed starting directly from the definition of 
the indicator function of the reaction event ^'at in Eq.(^. To do this, suppose first that for 
a given realization of traps' placement the nearest to the origin trap B appears at a certain 
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distance L. For such a realization, evidently, \&Ar = 1 for those trajectories {rjv} which do 
not leave within first steps the volume L'^. Consequently, here the overall decay function 
can be bounded from below by 

Pn > maxL[ exp [ - cL'^] Prob(|r„| < L for any n G [0, A^]|0)}, (19) 

where the first multiplier gives the probability of having a trap-free void of volume L'^, while 
the second one stands for the probability that a random walk starting at the origin does 
not leave this volume within first steps. For N sufficiently large, the latter probability 
follows exp — '-)clN/L?' , and hence, maximizing the rhs of Eg . ([191) with respect to L, i.e. 
searching for the maximal lower bound, one ends up with the dependence of essentially the 
same form as that given by Eg. ([18|) . 

Conseguently, these bounds demonstrate that the long-time behavior of the survival 
probability in Eg.(|^) is supported, first, by the presence of sufficiently large trap-free voids 
of typical size ~ A^i/{'^+2)^ second, by realizations of random walks constrained not to 
leave these voids within time A^, i.e. atypical, spatially-confined realizations for which |r7v| 
grows in proportion to A^^/('^+2) only. 

To conclude this section, we emphasize that the long-time behavior of the trapping and of 
the target annihilation problems are different. Whereas the integral effective rate constant 
for the target annihilation follows the behavior of the expected number of distinct sites 
visited, the trapping decay contains all the higher moments of this characteristic property 
and tends at very long times towards the asymptotic form in Eg. (p!6D , i.e. Q^^^ ~ {N/ cY/('^+'^) 
as N ^ oo. Thus the decay patterns differ, depending on which of the two species in Eg.(|I|) is 
the mobile one. One has conseguently a counter example to the view that in reaction kinetics 
only the relative motion of the species, but not the individual movements are important. In 
what follows we intend to analyse how the reaction kinetics will differ depending on which 
of the species precisely is being gated. 
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III. MODEL I: AN IMMOBILE TARGET A AND RANDOMLY MOVING 
STOCHASTICALLY GATED SCAVENGERS B. 



We start our analysis of stochastically-gated pseudo-first-order reactions considering first 
survival of an immobile target in the presence of a concentration of mobile gated scavengers 
B - the model which again admits an exact solution. We suppose here that each gate can 
be in either of two states, one - active and the other one - blocked. In the active state 
the B particles are reactive, whereas the blocked or inactive state inhibits the reaction. 
The gate on any B particle is supposed to update its state at each moment of time, at 
random, and independently of the gates imposed on the other particles. The A particle gets 
"annihilated" at the very moment when any of the B particles visits it for the first time 
being in the reactive state. On contrary, if any of the -Bs visits the A particle being in the 
blocked state, both particles can harmlessly coexist with each other. 

More precisely, we specify the reactive ability of the fc-th scavenger, (A; = 1, . . . , K), at 
the n-th step, n = 0,1, . . . , N, by assigning to each of the B particles a random variable 
r]!^\ This random variable may assume two values - with probability p, in which case the 
scavenger is neutral with respect to the reaction, and the value 1 with probability l—p, which 
corresponds to the reactive state. In all the models to be studied here, we will suppose that 
the reactivities of gated particles follow independent Poisson processes ||I4| , such that all rj^^ 
are independent, delta-correlated with respect to n and k randomly distributed variables. 
The average with respect to the distribution of rj!^^ will be denoted by the overbar. 

Further on, denoting as S^'^-' the position at which the k-th. scavenger appears at the n-th 
step for a given realization of its A^-step trajectory {S^^}, we construct the reaction event 
indicator function \E'iv, which now takes the form 



"^N = lim exp 



N K 



n=Ofc=l 

= n n exp [-/? J (S(f )) r^i^')] , (20) 

ra=0 k=l 

where is the indicator function showing whether the fc-th scavenger is at the lattice 
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origin at the n-th step, or elsewhere, Eq.(^. 

The indicator function \I/7v is equal to one if within the A^-step random walk none of K 
scavengers has visited the origin being in the reactive state, and turns to zero otherwise, i.e. 
in case when within the interval [0, A^] at least one of the scavengers has once visited the 
origin being in the reactive state. Its average over the reactivity fluctuations, i.e. the states 
of the gates 77^^^ can be performed very directly, since for the Poisson gating model under 
study all terms in the double product in Eq.(^) are statistically independent of each other. 
Consequently, we have that 

N K 



^ = n n exp [ - /5 T (si'^) 

n=0 fc=l 
N K . 

= n n (i-p)iimexp[-/5j(sw 

n=Ofc=l ' °° 



+ P 



(21) 



Further on, noticing that 



lim exp 



/?t(s«)] =i-j(s 



(fc) 



(22) 



and hence, that 



[1 — p) lim exp 



-/3J(S(^)) 



+ p = exp 



-ln{p), 



(23) 



we find that the indicator function of the reaction event "^n, averaged over the reactivity 
fluctuations, takes the form 



^N = exp[-apAf{{SP}) 



(24) 



where 



N K 

A/'({sS^}) = EE^(si'^^ 

n=0 A;=l 



(25) 



Note now that the functional M ({S^^}^ determines the number of times that K given A^- 
step random walk trajectories {S^"*}, with their starting points at fixed positions Sq^\ pass 
through the origin. Note also that here all walks contribute independently, which means 
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that the simultaneous visit of k walkers is counted k times, and hence, A/'(^{S^''}^ can be 
interpreted as the "residence time" at the origin for K independent random walkers. 

As a matter of fact, the number A/'(^{S^'*}^ can be also viewed from a different per- 
spective, which turns to be rather useful for illustration of the distinction between different 
models. Namely, on Fig. 4 we depict, for a one-dimensional case, a given realization of several 
random walks trajectories in form of "directed polymers" in 1 + 1-dimensions, the F-axis 
being the time n; in this language A/" ({Sj^''}^ appears as the total number of times that an 
n-axis is intersected on a segment [0, A^] by a brush of K phantom directed polymers with 
their ends fixed on a (i- dimensional substrate. The survival probability P/v can be thought 
off then as the generating function of the moments of the number of visits rendered by K 
independent walkers to the origin during "time" A^, or of the number of intersections of the 
n-axis on the segment [0, A^] by a brush of directed polymers. In the infinite space limit 
Pn is hence the generating function of moments of the number of visits to the origin by 
random walkers which are initially uniformly distributed with a given number density c on 
an infinite (i-dimensional lattice. 

We pass next to averaging over the scavengers' trajectories. Noticing that all multipliers 
in Eq.(|25|) with the same index k are again statistically independent of the multipliers with 
a different k, we may write down this average as 



k=l n=0 
K N 



k=l 



n=0 



1 f r ^ 

T7 E Eso <^ exp - ftp E J (S, 

So I n=0 



K 



(26) 



Turning to the infinite space limit, we find from Eq.(p6D that the survival probability P/v 
follows 



TV 



K 



E 



So i exp 



a 



n=0 



exp 



-c5:Es4M{o}({Stv})}], 



(27) 



16 



where A^{o}(^{SAr}^ is a functional of a given A^-step random walk trajectory, which is 
defined by 



N 



MM^n}) = [l-exp -«p5]X(S„) (28) 



n=0 



One notices now that A^{o}({SAf}j counts again a number of visits to the origin by a 
given A^-step random walk trajectory with its starting point being fixed at Sq. Namely, 
7Vl{o} (^{Sat}^ = for such A^-step trajectories which never pass through the origin, equals 
1 — p for such trajectories which visit the origin only once, and 1 — p' for those trajectories 
which do it exactly j times. Consequently, the expected value of the functional M.{o} ({Sat} 
which enters the exponential in Eq.(^), can be written down as the following polynomial 
in powers of the gating probability p: 

N 

Es4a1{o}({S^})} = E/5j^^(0|So)(l-p''^'^-°), (29) 

j=0 

where 5so,o is the Kroneker- delta and /5^^(0|So) is the probability that a simple random walk 
starting from the site Sq visits the origin in the first N steps exactly j times We adhere 
here to the definition for /3jj''(0|So) presented in Ref. p5|] and adopt the convention that the 
initial moment n = is regarded as the zeroth visit to the site Sq. 

We turn next to the analysis of the A^-dependence of the integral effective reaction rate 
for the model I, Q^^\ It follows from Eq.(p9D that Q^^^ is defined as the polynomial of the 
form 



QSi^ = EE/^J;?Ho|So)(i-j>^-+^^o.o) (30) 

So j=0 

To compute Q^^^ for any A^, we introduce the generating function 



oo 



Q^'\0 = E Qn^ = (1 -p) E A°^(o|o)e"^ + 

N=0 N=0 



CO 



+ T.[^-p'nT.P\i\me+ E E(l-p^jE/5)^\o|So)^^ (31) 

i=l N=l So,So^Oj=l N=l 

where we have made use of an evident fact that all /5^^(0|So) vanish for j > N. 
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Now, to evaluate the generating function Q^^\^) explicitly, we calculate three different 
sums entering Eg . (|3l|) . First of all, we find that 



oo oo A'' 

;i -p) E A°Ho|o)e^ = (1 - p){i + E (1 - E Fniome} 



N=0 



N=l 



n=l 
P. 



C CO oo -| 

:i - + - E Fn{o\o) y: = T^(i - Fm o) 

^ n=l N=n ^ ^ 



(32) 



^ n=l N=n ^ ^ 

where -F(0|0; ^) is the generating function of F„(0|0) - the average probability that a random 
walk trajectory starting at the origin returns to the origin for the first time exactly on the 



n-th step [^. Further on, we have 



j=l N=l 



-| oo 

T^T (1 - i^(0|0; 0) E(l - ^■^')^(0|0; 
J- ^ j=i 

l-pF(0|0;O(l+p(l-i^(0|0;O)) 
1-^ l_pF(0|0;O 



(33) 



and eventually. 



EE(i-p^)E/3)^^(o|So)e 

So7^0i=l N=l 

l-F(0|0;O 1 



N 



F{0\0;O 
1-e 



E(l-P^>(0|0;0^"' 



T.F{0\So;0 



P 



1-e (l-F(0|0;O)(l-pi^(0|0;O) 
l~p l-F(0|0;O r 1 



Fmo) 



P{0\0;O 



P(0|0;O 



l_^l_pF(0|0;OLl-e 
where P(0|0;^), as usual, stands for the generating function of the probability P/v(0|0) of 
having a walker at the origin on the A^-th step, provided that the walker started his random 
walk at the origin. Consequently, summing up the results in Eqs. (^) to (|^) we obtain the 
following exact expression for the generating function: 

(1-p) l-F(0|0;O 1 (1-p) 



(34) 



pF(o|0; (1 - V + (1 - p)Pm 

p 1 -1 



1 + 



(35) 



[i-p)pmo- 

S{^) being the generating function of the expected number of distinct sites visited by an 
A^-step random walk |^, S{^) = J2'n=o -^o{'5({SAr})|^^. In principle, the integral effective 
reaction rate for model I valid for any can be obtained now by inverting the result in 
Eq ■ (p5D , which will require, however, computation of very complex integrals. 



We turn now to the analysis of the large- iV behavior of the reaction rate . However, 
before doing it, it may be expedient to make first the following observation, which might 
seem to be quite surprising at the first glance: 

Namely, the result in Eq. (|35|) reveals that for recurrent random walks the leading large- iV 
behavior of the integral reaction rate Q^!^^ should be independent of the gating probability! 
It happens actually because for recurrent walks P(0|0; ^) — > +00 as ^ ^ 1^ {N ^ 00), and 
hence, the expression on the rhs of Eq.(P3D appears to be independent of p. In virtue of the 
Tauberian theorem, it implies that the effective integral reaction rate Q^^^ is independent of p 
when N 00. Moreover, the result in Eq. (|35D shows that for recurrent walks the leading at 
N ^ 00 behavior of Qtv^ is defined exactly by the expected number of distinct sites visited 
by an A^-step random walk. Consequently, for recurrent walks and N ^ 00 the target 
survival probability P/v is not influenced by fluctuating gates imposed on the scavengers and 
has exactly the same form for reactions which are subject to stochastic gating or reactions 
in which the scavengers are always in the reactive state [pO|-p2[|. 

On the other hand, such a behavior is not counterintuitive and agrees with our previous 
knowledge of the diffusion-limited reactions kinetics. The point is that imposing a fluctuating 
gate on otherwise perfect scavengers is in a way similar to imposing the constraint that 
annihilation of the target by a scavenger may happen with some finite probability, or at a 
finite rate prescribed by certain elementary reaction act constant K^i {K^i < 00). Following 
the seminal mean-field analysis of Collins and Kimball |3l| (see also Refs. and |^ for 



more details), the overall reaction constant taking into account both the constraints imposed 
by the elementary reaction act (finite K^i) and the transport limitations (in order to react, 
particles have first to find each other in the course of their random motions) follows 

{dQ^/dNy' = ^ + -^, (36) 

where Ksmoi is the so-called Smoluchowski constant which equals the diffusive current 
through the surface of an immobile, perfectly adsorbing sphere. Now, it is well-known 



[see e.g. Ref. ||18[ for more discussion) that for low dimensional systems the Smoluchowski 
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constant is not a real constant but rather a time-dependent coefficient which vanishes as time 
evolves. It means that in low dimensions random transport of particles offers progressively 
higher resistance with respect to the overall reaction rate than the constraints imposed by 
the elementary reaction rate, which results ultimately in the kinetics which is totally con- 
trolled by random transport of particles towards each other and is independent of K(.i. This 
is precisely the effect which we observe in case of low- dimensional stochastically-gated target 
annihilation problem. 

We note also parenthetically that similar effect has been predicted recently for low dimen- 
sional catalytically-activated binary reactions, in which case the particles' reactivity does 
not fluctuate in time but is rather a random function of the space variables [^]. It has been 
shown here that long-time kinetics is also insensitive to the concentration of the catalytic 
sites which promote reactions between randomly moving A particles and is independent of 
Kei- Of course, in higher dimensional space (such that d is greater than the fractal dimension 
of the random walk) the effective reaction rate does depend on the density of catalytic sites 
and Kf.1. Similarly, for stochastically-gated target annihilation reactions -P(0|0; 1~) is well 
defined for d > 2, which implies, by virtue of Eq.(P5D, that the leading at ^ oo terms in 
the integral effective reaction rate should depend on the gating probability p. 

We focus next on the special case of Polya random walks and proceed to determine 
the long-time behavior of the reaction rate in Eq.(^) explicitly, first for one- and two- 
dimensional lattices, in which case the Polya walks are recurrent, and then for (i-dimensional 
lattices with d > 2, for which the walks are non- recurrent. 



A. Polya walks on one-dimensional lattices. 

For Polya random walks on one-dimensional lattices the generating function of the first- 
visit probability is known exactly and has a particularly simple form (see e.g. Ref. [p5|]), 
F(0|0;^) = 1 — y/l — Consequently, in this case the generating function Q^^\^) of the 
integral effective reaction rate obeys 
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In the asymptotical limit ^ — > 1^ (or equivalently, when N — > oo) we find then from Eq.(p7D 
that 

'3'''«)-(r^-^(wfe) + o(Vv^)- (38) 

Hence, by virtue of a Tauberian theorem, we have that in the limit N ^ oo the effective 
reaction rate follows 



i.e., as we have already remarked, the leading behavior as — oo in case of the target 
annihilation problem with stochastic gates imposed on the scavengers is exactly the same as 
in case of its ungated counterpart, Eq.(|]). The first correction term, however, does depend 
on the gating probability p and diverges when p — > 1, i.e. in the limit when scavengers 
are being completely inert with respect to the reaction. Simple comparison of the first two 
terms in Eq.(^) shows that the universal, p- independent behavior is established when 
exceeds certain cross-over value A^*, such that A^* ~ 7rp^/2(l —pY- Note also that similar 



behavior has been predicted earlier in Ref. [0 within the framework of a continuous-space 
description. 



B. Polya walks on two-dimensional lattices. 

The generating function P(0,0|^) is not known explicitly for Polya walks on two- 
dimensional lattices. However, its asymptotical behavior as ^ — 1^ (or, equivalently, when 
A^ — s> oo) is well documented (see e.g. [^) and is given by 

p{m = ^ In (y^) (1 + o{i - 0), (40) 

where the constant C2 has been defined in the text after Eq.(|^), while the constant K equals 
4, 8 and 12 for hexagonal, square and triangular lattices, respectively. 
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From the latter equation we find than that the leading asymptotical behavior of Q^^\^) 
as — i> 1^ follows 



irC- 



2 



(l-02 1n(l-0 

_ , {l-p)HK) . CI . 1 ^ 

Hence, by applying the Tauberian theorem we find that for two-dimensional target annihi- 
lation with stochastically-gated scavengers the integral effective reaction rate obeys 

where 7 denotes the Euler constant. Note that again, in accord with our earlier prediction, 
the leading large- iV behavior appears to be independent of the gating probability and pro- 
ceeds exactly in the same way as for the ungated target problem. This long-time regime can 
be observed, however, at considerably longer times than that for the one-dimensional sys- 



tems; on comparing the first two terms on the rhs of Eq.(42) we infer that the corresponding 
cross-over time A^* is given by 

N* ^ exp f^^l , (43) 
II — pi 

i.e. is exponentially large when p — 1, while in one- dimensional systems this dependence is 
only algebraic. 

C. Polya walks on d-dimesional lattices, d > 2 

Lastly, we turn to the case of recurrent Polya walks, which case is realized, namely, for 
lattices with spatial dimension d > 2. Here the probability R of eventual return to the origin 
is finite, and consequently, we find from Eg . (|35|) that 



^^'^^^^ {i-p)pmi-)+p (i-o^^^Qi-o^/^^' ^^^^ 

which yields, in the large- iV limit, the following result 
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(,) ^ l-p 

(l-p)P(0|0;l-)+p 



(45) 



Hence, for lattices with d > 2 the decay of the survival probability is purely exponential in all 
dimensions. Note also that the exact result in Eq. (|4^) confirms in a way the mean-field result 
by Collins and Kimball, Eg . (|36D ; as a matter of fact, it appears that Eg . (145|) can be cast 
exactly into the form of Eg. (^) if we set K^i = {I— p) /p a.nd Ksmoi = 1/-P(0|0; 1~). Note also 
that our Eg. (|45|) confirms the conclusion of Szabo et al |ll|] concerning the possibility of the 
calculation of the steady-state stochastically-gated rate constant in terms of an appropriately 
defined ungated model. 



IV. MODEL II: AN IMMOBILE FLUCTUATING TARGET AND RANDOMLY 

MOVING UNGATED SCAVENGERS. 



We turn next to the survival of a stochastically-gated, immobile A particle, - a target, 
in the presence of ungated scavengers B, which perform independent random walks on a 
(i-dimensional lattice. For this model the indicator function of the reaction event can be 
written down as follows 



N 



lim exp 



N 



K 



(46) 



71=0 k=l 

where ?7„ is the indicator variable of the gate imposed on the target, while S^'^^ defines the 
lattice positions of the k-th. scavenger at the ra-th step, n = 0, 1, . . . , A^. We again suppose 
that the target reactivity assumes at random two values - 1 and with probabilities 1 — p 
and p, respectively. In the state ?7„ = 1 the target is accessible for reaction and can be 
annihilated by any of the scavengers arriving at the origin, while in the state ?7„ = reaction 
can not take place. 

Averaging first ^jv in Eg.(HB]) with respect to the fiuctuations of the reactivity, we get 



N 



= W hm exp 



n=0 

N 



fc=l 



N r r K 

n limexp -/3Ej(s(f)) 

n,=0 \ L k=l 



+ P 



(47) 
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Furter on, noticing that 



( 1 — p) lim exp 



k=l 



+ P 



1, Ef=iX(s(f))=o, 



(48) 



and hence, rewriting this expression as 



[1 — p) hm exp 



k=l 



+ p = exp 



-a,(i-jfx:j(s 



(fc) 



(49) 



we find that the indicator function of the reaction event, averaged over the fiuctuations of 
the target reactivity, attains the form 

(50) 



*A, = exp[-apAr* ({S;;)} 
where M* ({Si^'^}) is given by 

N / K \ 

A^*({s!J^}) = E(i-^ E^(s^lfO ) (51) 

n=0 \fc=l / 

Note now that the functional {l — X {^^^-^X {S^^^^^ ^ measures the occupancy of the origin 
at time moment n. It equals zero if none of K walkers is present at the origin at the time 
moment n and equals one if one or several scavengers appear at the origin at the n-th step. 
In this regard, A/'*(^{Sj^'*}^ is similar to the earlier defined functional H (^{S^"*}^ appearing 
in the analysis of the model I. Important difference is, however, that a simultaneous visit 
of the origin by several walkers is counted as a single visit, and consequently, M* ({S^"*}^ 
describes collective behavior of all K walkers, which can not be factorized, as it appears in 
the model I. This substantial distinction between the models involving ungated and gated 
targets has been noticed already by Szabo et al in Ref. 0| (see also Ref. for more 
details), who stated that the crucial difference between the case when the gates are imposed 
on i?s or on the A particle is that in the latter case the " switching of the A from the reactive 
conformation to a non-reactive one is felt simulateneously by all scavengers". It means, 
in particular, that if we define the functional M* ({Sj^^}^ using the "directed" polymers 
representation in Fig. 4, then it would count all sites on the n-axis visited simultaneously by 
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two, three and etc walkers as singly visited sites. In this regard, Af* ({Sy^^}^ determines the 
number of distinct visits to the origin by K independent walkers. 

We notice next that an upper bound on the integral effective reaction rate Q^^'' for the 
model II can be found very straightforwardly. To do this, it suffices merely to observe that 

\k=l ) k=\ 

and hence, that M* (^{Sj^^^ < M ({S^''}^ . This implies, in turn, that the survival probability 
Pat for the model II is greater than the survival probability obtained for the model I, and 
the integral effective reaction rate Q%^^ obeys 

QSi'^ < (53) 

which inequality sets a rigorous upper bound on Q%^^ . Note that the inequality in Eq.(p3D 



has been established earlier using different type of arguments in Ref. p5| . 

It may be worthy to remark that the inequality in Eq.(|52D holds as an equality when 
all S^'^-' are different at a given n, which happens, namely, when the scavengers do obey a 
hard-core exclusion and none two scavengers can simultaneously occupy the same lattice 
site. Given that the scavegers are indistinguishable, we may thus expect that for the model 
II with hard-core scavengers the A particle survival probability will be determined exactly 
by Eqs.(^) to ( ^5l) at sufficiently large times. 

We proceed next to calculation of the upper bound on the survival probability for the 
model II, which requires a bit more complicated analysis. To do this, we return to the 
indicator function of the reaction event in Eg . (^7]) and perform first averaging with respect 
to the trajectories and the initial positions of the scavengers. The steps involved in the 
averaging procedure in this case are as follows: 

K N K N 

<^^> = < n limexp{-/3Er^„j(s(f))}>= J] < jim exp { - /3 J] j(s(f)) } >= 

k=l'^^°° n=0 k=l n=0 

K -I N 

= II{mJ: Eso{ ji^m exp [- (3 Y^Vn x(s(^))] }) = 

k=l So n=0 

= (^EP^ob(r/„T(S„) =Ofor anynG [0,iV]|So))'', (54) 

So 
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where Froh(rin 2^(S„^ = for any n G [0, A^]|So) is the probabihty that an A^-step random 
walk, starting at Sq and characterized by internal two-state variable ?7„, has never visited 
the origin being in the reactive state ?7„ = 1. Next, turning to the infinite-space limit, we 
find 



M . 

»0 



K 



< > = (l - ^ E (l - Prob(r7„ J(S„) = for any n G [0, iV] |S 
= exp 

Evidently, 

(l - Prob(r/„ t(S„) = for any n G [0, iV]|So)) = Prob( ^ Vn x(S„) > l|So), (56) 



c E (l - Prob(77„ J(S„) = for any n G [0, iV]|So))J (55) 

So 



N 



n=0 



where PTohyJ2n=oVn 2^(S„j > l|Soj is the probability that an A^-step random walk starting 
at the site So has at least once visited the origin being in the reactive state. 

Note that PToh(^J2n=oVn 2^(S„) > l|So) is not constrained in the sense that it provides 
no information at which of the visits to the origin precisely the reactive state has appeared; 
that is, the particle could visit the origin many times until it arrived eventually being in the 
reactive state. Having this in mind, we now realize how the sum in the exponent in the last 
line of Eq.(|55D can be bounded from below, which will result in the desired upper bound 
on the target survival probability in the model II. To do this, we will proceed as follows: 
Suppose first that the starting point of the trajectory {Sjv} is not the origin. Then, we 
notice that for So 7^ 0, one has 



N N / n-1 \ 

Prob( E Vn X(S„) > l|So) > E Prob rin x(S„) = 1 and ^ Vi ^(Sz) = 0|So 

n=0 n=0 \ 1=0 I 



(57) 



where Prob (^r^^ X^^S^^ = 1 ?i.vAYi!iZl Vi'^i^i) = 0|So) stands for the probability that the 
origin has been visited for the first time exactly at the n-th step (has not been visited prior 
to the n-th step) and the particle at this very step was in the reactive state, i.e. such that 
rjn = 1. Hence, the rhs of Eq.(|57|) is the constrained probability that within an A^-step walk 
starting at the site So 7^ the particle happened to be in the reactive state at its first visit 
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to the origin. Summing next both sides of the inequahty in Eq. (|57|) over all starting points, 
we obtain 

^Prob(^r7„ J(S„) > l|So) > E Prob( r^^ j(S„) > l|So) > 

So n=0 So.SoT^O n=0 

N / n-1 \ N 

> EP^ob h„T(S„) =land5]r/, J(S,) =0|S 

So,So7^0n=0 \ 1=0 

N N N 

n=0 So.SoT^O 



E E^nFn(0|So) 

So,So^On=0 



(58) 

n=0 So n=0 

where A(^{S„}^ is an auxiliary random variable, defined to be the number of virgin sites 
visited on the n-th step by a given particle trajectory {S„} p5| , 

A({SJ) = 5({SJ) - 5({S„_i}) (59) 

Consequently, we can bound the rhs of Eq. (F^) as 



Pn < exp 

N 



c E^n^o{A({S4)} = nexp[-cr/„Eo{A({Sj)} 

n=0 J n=0 

n {(1 -p)exp[- cEo{a({S4)}] +p] 



(60) 



n=0 



> . . . > 



i?o{A({S^})}, (61) 



(62) 



Further on, following the Dvoretzky-Erdos Lemma -Eo{a(^{S„}^ | is a monotonic decreasing 
sequence of time n (see, e.g. [^), i.e. 

1 > Eo{a({Si})} > i?o{A({S2})} > . . . > i?o{A({Sj) 

such that 

/^m™i?o{A({S4)} = -^^^^ 

we can majorize the terms in the curly brackets on the rhs of Eq.(^0]) by replacing 
-E'o{a(^{S„}^ I by its minimal value i?o|A(^{SAr}^ |. Enhancing in such a way the bound in 
Eq.(|60D, we have that 

N 

Pn<X{ {(l-p)exp[- cEo{a({S^})}] +p] = 

n=0 

'l-p)exp[-cEo{A{{SN})}]+py 



N+l 



exp 



{N + 1) In (l/((l - p) exp[- c Eo{a({S;v}) }] + p)) 
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(63) 



d=2 gSi^)> (l-p)f^(l + 0(l/ln(iV) 



and hence, the integral effective rate constant obeys 

Q(i^)>iil±l) ln(l/((l-p)exp[-ci?o{A({SA.})}]+p)) (64) 

For (i-dimesnional Polya walks, in particular, we find from Eq.(]6^) the following explicit 
lower bounds on the integral effective rate constant in the model II: 

d=i Q'i''> {i-p){^y''{i+o{i/VN' 

In(A^) 

d>2 Qf^>^ ln(l/((l-p)exp[-c/P(0|0;l-)]+p))(l + O(l/yiV)), (65) 

which hold in the limit N ^ oo. 

On comparing the results in Eqs. (|53|) and (|65|), we notice that both lower and upper 
bounds display the same iV- dependence, but differ slightly in numerical factors. This means, 
in turn, that the iV-dependence of the integral effective reaction rate Q^^^ is determined 
here exactly by Eqs. (|53D and ([651). Consequently, at sufficiently large times the decay laws 
in the models I and II of gated target annihilation are essentially the same (up to a possible 
difference in characteristic decay times), and coincide with the decay law predicted for the 
ungated model of Section II. 

Note also that the time evolution of the A particle survival probability for the model 
II defined in a one-dimensional continuum has been considered earlier in Ref. |]13|. Within 
the framework of the heuristic Smoluchowski-type approach, it has been predicted that the 
long-time decay of should follow the decay law in the first line in Eq.(H), i.e. should 
proceed at long times exactly in the same fashion as that for the model I and consequently, 
should be essentially the same as in the ungated target problem. While intuitively such a 
behavior seems to be quite plausible for low- dimensional systems (see the discussion following 
the Collins-Kimball result in Eq.(p6D) and, as a matter of fact, agrees with our prediction 
in Eq.(^), one still can not, in principle, rule out the possibility that the integral effective 
reaction rates for models I and II may have different numerical factors even in low dimensions. 



The point is that the Smoluchowski approach in Ref. |T^, which is a continuous-space 
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counterpart of the Rosenstock approximation, allows to determine here only a lower bound 
on the target survival probability and thus can not produce exact numerical factors. 



V. MODEL III: AN UNGATED, MOBILE A PARTICLE AND IMMOBILE, 
RANDOMLY PLACED, STOCHASTICALLY-GATED TRAPS B. 



We turn next to the case of stochastically-gated trapping reactions, focusing first on 
the situation involving an ungated A particle, which performs a discrete-time, homogeneous 
random walk on a d-dimensional lattice starting from the origin at n = 0, in the presence 
of immobile, randomly placed, stochastically-gated B particles. The properties of the gates 
are the same as defined in the model I. 

For this model, the indicator function of the reaction event can be written down as 
follows 

N K 



= lim exp 



(fc) 



N K 

— Y\.Y\. ^™ 



n=0 fe=l 



(66) 



n=0 fe=l 

where 8*^*^^ are d-dimensional lattice vectors denoting positions of K traps B, while r„ defines 
the lattice position of the A particle at the n-th step. 

Averaging first over the reactivity of the traps, we readily find 



N K 



^iv = n n exp -(31{vn - S(^')) n^^^ 

N K 



n=0 k=l 

Further on, noticing that 



n n {(1 - p) j™,e^p [ - /3 2: (r, - s^'^)) ] + p] 



(67) 



{1 — p) lim exp 

/3— »oo 



-/3X(r„-S«) 
and hence, rewriting this expression as 



+ p 



r = S(*^) 



[1 — p) lim exp 

/3— »oo 



-/?j(r,-S^' 



p = exp 



P 



1, r„^SW, 



- a^j(r, - S^) 



(68) 



(69) 
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we have that the indicator function of the reaction event, averaged over the reactivity fluc- 
tuations, takes the following form: 

N K 



"^N = exp 



(70) 



n=0 fc=l 

Note, that the averaged indicator function ^'at is now an exponential of the factor ap times 
the number of times a given random walk trajectory starting at n = at the origin visits a 
given array of lattice sites {S'-'^^}, i.e. can be thought off as the moment generating function 
of the "residence time" of a single random walker on a given array {S'-'^-*}. From a different 
perspective, can be viewed as the partition function of a phantom polymer chain on a 
lattice with randomly placed energetic barriers of finite height; the limit ap — > oo, {p —>■ 0), 



corresponding then to the case of randomly placed, impenetrable obstacles 



Now, the double average over the trajectories of the A particle and over the positions of 
the traps can be written down as follows 



Pn=< > =Eo{< nexp[-ap ^ J (r„ - S^'^)) ] > } 

k=l n=0 
K N 

= Eo{n <exp[-ap ^J(r„-S(^'))] >} = 



fc=i 



n=0 

N 

ap J2 ^{^n - S 



K 



} 



(71) 



S n=0 

Turning next to the infinite-space limit, i.e. letting M, K oo, while keeping their ratio 
fixed, we have that 



1 ^ 
P^ = i?o{(l-^E(l-exp[-«p E^( 



r„ - S 



n=0 



= Eo{exp -c^A^{s}({rjv}) }, 
L s J 

where A^{s}({rAr}) is defined by 

/ N 
A<{s}(r„) = 1 - exp -ap^J(r„-S) 



(72) 



(73) 



n=0 



Note that 7Vl{s} ({r7v}j is quite similar to the functional defined in Eq.(|2^) with the only 
minor difference being that the latter is associated with the multiple visits to some given site 
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S; namely, A1{s}({rjv}^ = 1 — if the site S is visited exactly j times by an A^-step walk 
starting at the site tq = 0. Consequently, the sum Z^s -^{s}({i"Ar}) probes the occupancy 
of the sites visited by a given realization of a random walk trajectory (see Fig. 5). 

However, an important difference with the previously considered model, which makes 
the computation of P/v for the model III to be a fairly complex problem, is that here we 
have to deal with a moment generating function of A^{s} ({rAf}) , compared to a much easier 
problem of computation of an expected value of 7W{s}(^{r7v}), encountered in the model I. 
Below we consider first predictions based on some approximate approach - the Rosenstock 
approximation, and next derive rigorous lower and upper bounds, which both have the same 
dependence on the time but differ in the prefactors. 

A. Rosenstock approximation. Decay pattern at intermediate times N. 

We start our analysis of the A^-dependence of the survival probability in Eg . (|72D by 
considering first the predictions of the Rosenstock approximation. Applying the Jensen 
inequality, we may bound P/v in Eq. (|72|) as follows 

Pn > exp [ - cEo{Y.M{s}{{rN})}] (74) 



To proceed further, we note that the sum I]s -^{s}(^{rAr}j can be written down formally 
(see Fig. 5) as the following polynomial with random coefficients 

N N 

E^{s}({rA.}) = E (I-P^) = 5({r^.}) - E V^^'^ ({r^}) ^ (75) 

s j=l j=l 



where each V^^^\^{rN}j, j = 1, N, is a random variable, which equals the number of sites 
in a given iV-step random walk trajectory {r^} being visited exactly j times. Note that 
V*^-'-* (^{ttv}) are not independent; this can be readily seen if one notices that the combination 
SjLi jV*--^''(^{rjv}^ is a non-fiuctuating quantity and equals the total number of sites visited 
by an A^-step walk, i.e. A^. Lastly, the random function iS(^{riv}) in Eq.([73D denotes, as 
previously defined, the number of distinct sites visited by an A^-step random walk trajectory 

{M,5({M)=Ef=iV(^)({M). 
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Consequently, we find that Eg. ([7^) attains the following form 



N 



Pn > exp [ - cEo{5({rjv})} + Eo{v(^) ({r^v}) } p^] (76) 

The averaging in the exponent in Eq.([76|) can be performed directly using the results ob- 
tained in ||2^ for the generating function of the expectation Eq^V^^^ (^{rN}^^ the other 
hand, it is evident that for homogeneous random walks 

^o{E>^{s}(W})} = E^r„{A^{0}({rAr})}, (77) 

S ro 

which implies, in turn, that in terms of the Rosenstock approximation the integral effective 
reaction rate for the model 111, Q^l/^\ has exactly the same form as that obtained for the 
model 1. As a matter of fact, this does not seem to be surprising since the mean-field- tj^e 
Rosenstock approximation is not sensitive to the fact which of the reactive species precisely 
is mobile and which is fixed (see Eqs. ([7^) and (0)). However, a profound difference between 
these two models does exist and below we will show that the large- decay of P/v proceeds 
slower than it is predicted by Eqs. (39), (42) and (45). 



B. Large- decay of the survival probability. 



First of all, we note that in virtue of Eq. (|7^) and of an evident observation that for any 
given random walk trajectory {rjy} all V'--'^ (^{r^v}) are non-negative, one finds that 

-M{s}({r^}) <5({r^}), (78) 

which implies that, quite trivially, P/v is bounded from below by 



P^>Eo{exp[- c5({r,v})]} 



(79) 



On the other hand, for any p <1 and any j > 0, the difference l—p' > 1—p. Consequently, 
the survival probability is always bounded from above by 



Pn < Eo{ exp [ - c (1 - p) 5({r,v})] } 
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(80) 



i.e. large-A^ decay of Pat in the model III proceeds slower than the decay in the ungated 
trapping problem (section II) with the concentration of traps equal to c(l — p). This in- 
equality, however, does not seem to be trivial at the first glance, since for the model III the 
factor c(l — p) represents only the average value of the active traps concentration, which 
does fiuctuate in time and may exceed c(l — p) at certain time moments. 



Lastly, taking advantage of the analysis in Ref . |^ , we infer from two latter inequalities 
that the integral effective reaction rate for the model III obeys the following two-sided 
inequality: 

a, (I- rt^/<^«) {^f*'> < < a, {^f''*-" (81) 

Note now that both lower and upper bounds show the same dependence on A^, (but have 
slightly different prefactors), which insures that in the large- limit the integral effective 
reaction rate Q'^}/^^ grows in proportion to A^'^/('^+^) and consequently, the decay of the 
survival probability for the model III is described by the dependence ln(P/v) ~ — A^'^/('^+2)^ 
i.e. the same iV-dependence as in the ungated case 



VI. MODEL IV: A MOBILE, GATED A PARTICLE AND RANDOMLY PLACED 

FIXED TRAPS 

Consider lastly the trapping model involving an A particle bearing a stochastic gate and 

performing lattice random walk in the presence of randomly placed, immobile ungated traps. 

For such a model the indicator function of the reaction event can be written as 

[ N K n 




which function equals unity if an A^-step trajectory {rjv} does not visit any site from a given 
array {S'-'^^} being in the reactive state, and turns to zero if any of r„, n = 0, 1, . . . , A^, 
coincides with any of S*-'^-' when = 

We turn first to averaging the indicator function of the reaction event, Eq.(^2]), over the 
traps' placement on the lattice. Rewriting first \l/7v in Eg . (|8^) in the factorized form 
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K 



= W lim exp 



fc=i 



N 

I 

?i=0 



(83) 



and noticing that since all traps are placed independently of each other, all multipliers in 
Eq.(p3D appear to be statistically independent, we have that the indicator function averaged 
over the traps' placement reads 

N 



1 

*^){s«} = (m Eji^exp [ - (3 Y.nnl{r^ - S 



K 



(84) 

S n=0 

where the brackets with the subscript {S'-'^^} stand for the averaging with respect to the 
positions of the traps. Next, turning to the infinite-space limit, we find that (^^n^ is 
given explicitly by 



{SW} 



N 



{SW} 



exp{ - cS({YN}\{r]n = 1})}, 



(85) 



where the functional 

S{{rN}\{Vn = 1}) = E (l - ^( E Vnl{rn " S 



N 



(86) 



S n=0 

determines the number of distinct sites visited in the rjn = 1 state by a given A^-step 
trajectory {tn}, or, in other words, the number of distinct sites visited by the A particle 
being in the reactive state within a given realization of its A^-step random walk. Below we 
will study the temporal evolution of the function in Eq.(|85D, averaged over the reactivity 
fluctuations and trajectories {rjv}, using first, the Rosenstock approximation, and then, by 
evaluating rigorous lower and upper bounds. 



A. Rosenstock approximation. An intermediate time decay. 

Consider now the prediction of the Rosenstock-type approximation for the A particle 
survival probability in the model IV. Applying the Jensen inequality, we have then that the 
particle survival probability in the model IV is bounded by 



PN>exp[ - c Eo{5({rjv}|K = 1})} 
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(87) 



The average of S(^{T]\f}\{rin = 1}) over the reactivity fluctuations can be performed straight- 
forwardly. First of all, we rewrite Eq. (|86D as 



N 

S{{rN}\{Vn = 1}) = E (l - ^( E Vnl{rr. " S))) 

S n=0 

1 /-s^ 
27r io 



N 

dZ W exp iZr]nl{rn - S) ) 



(88) 

Next, averaging the latter equation with respect to the distribution of the variables {rjn}, 
we have that 



1 



S{{rN}\{Vn = 1}) = E (l - ^ X + (1 - 

^E-^s({M), 



iZ 



(89) 



where the functional A^s({rAr}) has been defined previously in Eq. ([75D . 

On comparing the result in Eq.(pUD with Eqs.(0) and (^), we notice that within the 
Rosenstock approximation the integral effective rate constant for the model IV appears to 
be exactly the same as the one previously found for the model III and coincides, as well, 
with the result obtained for the integral effective rate constant in the model I. Hence, this 
approximation predicts that the decay of the survival probability in the model 111 proceeds 
exactly in the same fashion as the decay laws obtained for the models 1 and 111. In other 
words, the Rosenstock approximation appears to be completely insensitive to the fact which 
of the reactive species precisely is mobile and which of them precisely is being stochastically 
gated. 



B. Large- decay of the survival probability. 

We start with the derivation of a rigorous lower bound on the A particle survival proba- 
bility. Here, averaging of the indicator function of the reaction event in Eg . (|8^) with respect 
to the reactive state of the mobile A particle can be performed as follows: 



N 



lim exp 



TV K 
-(3 E^n E^(rn-S«) 

n=0 k=l 
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TT lim exp 

n=0 

N 



K 



k=l 



N f r X 

n (1 - lim exp -PY.^ (rn - S^'^) 

"1=0 k L fc=l 



(90) 



Next, noticing that 



'1 — p) lim exp 



K 



fc=i 



1, Ef=iX(r„-SW) =0, 
p, otherwise, 



(91) 



and hence, that 



[1 — p) lim exp 



-/5X:x(r„-S 



(fc) 



fc=l 



+ p = exp 



«p(l-^fE^(rn-S(^))))], (92) 



we find that the averaged indicator function of the reaction event obeys 



where 



'^N = exp [-Op A/'{s(fc)}({rAf})], 



n=0 \ \fe=l / / 



(93) 



(94) 



is the "residence time" of a given A^-step random walk trajectory on the subset of "distinct", 
i.e. non-coinciding sites from the set {S*-^^}, k = 1, . . . , K. This means that if any two (three 
and etc) sites from {S*^'^^} coincide, i.e. the traps overlap, a visit of {r^v} to such a multiply 
covered site singly contributes to A/|s(fc)} (^{rjv}^ . A rigorous lower bound on \l>7v follows then 
from an evident inequality 



N / / K \\ N K 

E 1 E^(rn - s(^)) < E E I (r. - s^'^) 

n=0 \ \A;=1 // n=Ofc=l 



(95) 



where the rhs determines the unconstrained "residence time" of the same A^-step random 
walk trajectory on the set of all sites from {S'''^^}. Clearly, the inequality in Eq.(|95|) becomes 
an equality if all sites {S^'^-'} are distinct, e.g. when the traps do obey a hard-core exclusion. 
Consequently, in virtue of the inequality in Eq. (|95D , we have that 



N K 



^^>exp -ap^^j(r„-S(^) 



n=0 fc=l 



(96) 
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where the rhs, as one can readily notice, is exactly the indicator function of the reaction 
event for the model III, averaged over the reactivity fluctuations. This implies the following 
inequality 

or<«r <«.(^)'"""', (97) 

which signifies, in particular, that similarly to the relation between two Qn^, describing 
survival of gated and ungated targets (models II and I), the integral effective reaction rate 
for the model involving a mobile gated A particle in the presence of immobile ungated traps 
is generally less or equal to the corresponding rate for the model with ungated A particle 
performing random walk in presence of gated traps. 

We proceed finally to the derivation of a rigorous upper bound on the A particle survival 
probability in the model IV. To do this, let us turn back to the functional S(^{rN}\{rin = 1}) 
in Eq.(p6|) and recall that it determines the number of distinct sites visited by a particle 
appearing in the reactive state within its A^-step random walk {ttv}. Note now that similarly 
to the situation encountered in the derivation of the analogous bound in the model II, here 
we do not have any restriction at which visit precisely the particle appeared in the reactive 
state; that is, each site S can be visited many times by inactive particle until it re-appears 
at this site eventually being in the reactive state, which event only does contribute to 
the overall value of the functional S(^{r]\f}\{rin = 1}). Hence, to find a lower bound on 
S(^{rN}\{rin = 1}) in Eq. (|8^) we will pursue the strategy employed already in the Section 
IV, i.e. we will try to restrict the order of the reactive visit to the site S. Here, however, 
it appears to be a bit more delicate problem, since we have to deal with the realization- 
dependent functional in Eg . (pB]) , rather than with its expected value. 

To find a lower bound on S(^{rj\f}\{rin = 1}), we thus proceed as follows: First, 
we constrain the summation in Eq. (|86|) supposing that it runs not over all sites of an 
infinite lattice, but only over some finite subset {S*}. Clearly, since the functional 
(^1 — l(^J2n=oVnl{T'^n ~ '^))) positive-dcfinite, such an operation will result in a lower 
bound on s(^{rN}\{rin = 1}). Next, we define the subset {S*}; we stipulate that for a given 
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realization of particle's trajectory the subset {S*} is just a collection of such lattice sites S, 
on which the particle appeared for the first time being in the reactive state, i.e. sites which 
remained "virgin" until the first visit in the rjn = 1 state. 

More formally, derivation of such a lower bound on S(^{rN}\{rin = 1}) can be based on 
the evident inequality between the following two indicator functions: 

N 

Is{{rN}\{Vn = 1}) = (l - X( E Vnl{r^ - S))) > 

n=0 

N n—1 n 

> As({r^}|K = 1}) = E E ^(r^ " S)) - ^( E^(r/ " S))}, (98) 

n=0 1=0 1=0 

where the indicator function on the Ihs of Eq. (|98D equals one if the site S has been visited 
at least once by the particle being in the reactive state within a given realization {rjv} of its 
A^-step walk, and zero - otherwise; while As({rAr}|{?7„ = 1}) equals 1 if within an A^-step 



walk the particle appeared for the first time on the site S being in reactive state not having 
visited this site before; and equals zero otherwise. 

Summing both sides of the inequality in Eq.(p8D over all lattice sites S, we have, conse- 
quently, that 



S{{rN}\{Vn = 1}) = E^s({r^}|K = 1}) > 

s 

N n-1 n 

> E E^4^( - S)) -l{T.^{ri - S))} = 

S n=0 1=0 1=0 

N n n-1 

= E^^E {(l -^(E^(r. - S))) - (l -X( E^(r. - S)))}, (99) 

71=0 S Z=0 1=0 

which yields, in virtue of the definition in Eq. (|59|) , the desired lower bound of the form 

N 

5({r^}|K = 1}) > E^nA(K}), (100) 

n=0 



where the rhs of Eq. (|100|) determines the number of "virgin" sites visited by an A^-step 



random walk being in the reactive state. Equation ( |100|) implies, in turn, that the function 



^^/{sc^)} Eq.(|85|) is bounded from above by 



N 

^^>rs(.), < n exp [ - c A({r„})] (101) 

^ ' n=0 
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Now, averaging both sides of the latter equation over the fluctuations of the reactivity, we 
find 



N 



N N 

< n exp -c rin A({r„,}) = exp -c rjn A({r„}) 



n=0 



n=0 
N 



n {(1 - P) exp [ - c A({r„})] + p} (102) 



71=0 



Recollecting next that the realization-dependent property A(^{r„}j assumes only two values 
- 1 or 0, and hence, that 



l-p)exp -cA({r„}) +p 

"1, A(K}j=0, 
we may rewrite quite formally the bound in Eq.( |102| ) as 



{(l-p)exp[-c]+p}, A({r„}) = 1, 



(103) 



N 



{SC^)} 



N 

< n {(l-p)exp[-cA(K})] +p} 



71=0 

exp 



N 



In 



il-p)exp[-c]+pj ^ 



(104) 



Lastly, noticing that Z]^=o^({^n}) = '^^^jr^v}), we find that the A particle survival proba- 
bility obeys 



In 



1 



5({rAr} 



}, (105) 



^1 -p)exp[-c] +p^ 

and hence, by taking into account the lower bound in Eq.(p7|), we arrive at the following 
two-sided inequlity for the integral effective rate constant in the model IV: 

2/{d+2) 

^./(.+2) < QilV) < ( 

C 



ad 



In 



(106) 



'1 - p) exp[-c] +pj J 

Note that again, both sides of the inequality in Eq. ( |106| ) show the same dependence on the 
time N and consequently, determine exactly the A^-dependence of the integral effective rate 
constant in the model IV. We also remark that in the limit c <^ 1 one has that ln(l/ ((1 — 
p) exp[— c] +p)) ~ (1 —p)c and hence, in this limit the lower bound on the integral effective 
reaction rate for the model IV, Eq.( |106| ), coincides with the lower bound on Qtv^^"* in Eq. (pl|) . 
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C. Conclusion 



To conclude, we have studied the time evolution of the A particle survival probability in 
four models of stochastically-gated, diffusion-limited pseudo-first-order reactions of the form 
A + B ^ B. We have considered two different models of target-like annihilation reactions, 
where the first one concerns the survival of a single, immobile A particle (the target) in 
the presence of randomly moving gated scavengers B (model I), while the second model 
focuses on the fate of a gated immobile A particle in the presence of randomly moving 
ungated scavengers B (model II). Two other examples of stochastically gated pseudo- first- 
order reactions are furnished here by the trapping reactions between a mobile, ungated 
A particle and a concentration of randomly placed, immobile, gated traps B (model III) 
and the reverse situation with a mobile gated A particle and randomly placed, immobile, 
ungated traps B (model IV). In all the models studied we have supposed that the mobile 
species perform symmetric lattice random walks. Besides, we have adopted here the two- 



state Poisson gating model of Ref . |]T4| , in which each of the gates is supposed to be in either 
of two states - an active and a blocked one, and to update its state at each tick of the clock 
at random, independently of the previous history as well as of the gates imposed on other 
particles. 

We have demonstrated that the model I allows for an exact solution and derived the 
explicit asymptotic decay forms for lattices of different spatial dimensionality. Curiously 
enough, it appeared that for low dimensional lattices {d < 2) for which the lattice random 
walks are recurrent, the long time behavior is independent of the presence of stochastic gates 
(as soon as the gating probability p < 1) and proceeds exactly in the same fashion as for 
reactions with well-defined, non-fluctuating reaction rates (Section II). Correction terms do, 
however, depend on the gating probability p and may be important for reactions in which 
the species are being blocked most of the time. Next, we have found that for the model I in 
higher dimensions the decay is described by a purely exponential function of time with the 
characteristic time dependent on both the gating probability and on the probability of the 
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eventual return to the origin. Physical explanation of the predicted behavior has been also 
provided. 

Further on, for the model II the decay pattern has been determined exactly in form of 
rigorous lower and upper bounds showing the same dependence on the time N but having 
slightly different prefactors. We have demonstated that the decay of the A particle survival 
probability in this model is characterized by essentially the same time dependence as that 
for the model I, i.e. the integral effective reaction rate follows the behavior of the expected 
number of distinct sites visited by an N-step random walk but may have a different numerical 
factor. 

Next, for models III and IV we have presented some approximate results, based on 
the so-called Rosenstock approximation which may provide a plausible description of the 
kinetic behavior at intermediate times, as well as exact results concerning the long-time 
evolution of the A particle survival probability. We have demonstrated that within the 
Rosenstock approach no difference exists between the kinetic behavior in the models III and 
IV; Moreover, we have shown that the decay forms coincide with the exact result obtained 
for the model I. Long-time evolution of the decay functions in the models III and IV has 
been determined in form of rigorous lower and upper bounds characterized by the same 
dependence on the time A^. We have also realized that in the case of stochastically-gated 
trapping reactions the long-time decay of the A particle survival probability has essentially 
the same form as that describing the kinetic behavior of their ungated counterparts (Section 
II); the characteristic times might be, however, renormalized to include the dependence on 
the reaction probability. 

As an interesting by-product of our analysis, we have also shown that the survival prob- 
ability in all four models under study can be interpreted as a moment generating function 
of some refined characteristics of random walk trajectories. In particular, we have demon- 
strated that for the model I the survival probability is the moment generating function for 
the number of visits rendered by a concentration of independent random walkers to the ori- 
gin. In other models this survival probability appears as the moment generating function of 
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the number of self-intersections of random walk trajectories, residence time on a disordered 
array of marked sites, the number of sites visited exactly a given number of times and so on. 
Consequently, our results apply as well to the asymptotical behavior of the above mentioned 
generating functions, which in many cases has not been known yet. 
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Figure Captions. 



Fig.l. A schematic illustration of effective geometrical screening of an active A particle 
by inactive parts of a complex (polymer) molecule. In the situation depicted on this figure 
the A particle, which is attached to a polymer, is completely inaccessible to the B species 
and the reaction between them is inhibited due solely to the geometrical restrictions. 

Fig. 2. In the situation depicted in this figure the Bs may diffuse through the hole 
opened in the course of polymers' thermal motion and hence, may enter into reaction with 
the chemically active A particle. 

Fig. 3. Pseudo-first-order A + B ^ B reaction involving a single mobile (immobile) A 
particle and a concentration of fixed (mobile) B particles taking place on a two-dimensional 
lattice. The gate may be imposed on either A or B particles. 

Fig.4. Representation of the B particles trajectories in form of "directed polymers". 
Numbers on the n-axis denote the total number of crossings of different points on this axis 
by different trajectories. 

Fig. 5. A realization of N — 130-step random walk trajectory {rjv} on a two-dimensional 
square lattice. The sites visited two times are marked by circles; the sites visited 
three times - by squares and four times - by diamonds. For this particular realization 

A^{S}({riv}) is a fourth-order polynomial with respect to the gating probability p of the 
form At{s}({rAr}) = 113 - 69p - 28p^ - 12p^ - 4p^. The coefficients in this polynomial 
are random, correlated variables dependent on the particular realization of trajectory {rjv}. 



45 




Fig.l. Benichou et al., PRE 




Fig.2. Benichou et al., PRE 
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